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Electron Localization in a 2D System with Random Magnetic Flux 
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Using a finite-size scaling method, we calculate the localization properties of a disordered two- 
dimensional electron system in the presence of a random magnetic field. Below a critical energy E c 
all states are localized and the localization length £ diverges when the Fermi energy approaches the 
critical energy, i.e. £(E) oc \E — E c \~ lJ . We find that E c shifts with the strength of the disorder 
and the amplitude of the random magnetic field while the critical exponent (y ~ 4.8) remains 
unchanged indicating universality in this system. Implications on the experiment in half-filling 
fractional quantum Hall system are also discussed. 

PACS numbers: 72.10.Bg, 71.55.Jv, 72.15.Rn, 73.40.Hm 



Fractional quantum Hall (FQH) system is an ideal 
candidate to study localization properties in strongly cor- 
related electron systems. In a non-interacting picture, 
according to the scaling theory of localization [H all elec- 
trons in a two-dimensional system are localized in the ab- 
sence of magnetic field. When the two-dimensional elec- 
tron system is subject to a strong perpendicular magnetic 
field, the energy spectrum becomes a series of impurity 
broadened Landau levels. Extended state appear in the 
center of each Landau band, while states at other ener- 
gies are localized. This gives rise to the integer quantum 
Hall effect. It has been demonstrated that the critical 
transitions in the center of each Landau band are univer- 
sal ||. In the FQH regime, electron-electron interaction 
plays an important role Q|, however, recent experiments 
and theories indicate that the critical properties of the 
plateau transition might also be in the same universality 
class §,§. 

Recently, Halperin, Lee and Read j|] an d Kalmeyer 
and Zhang developed an effective Chern-Simons field 
theory to understand electronic properties of the FQH 
systems. In their theory the quasi-particlcs are weakly 
interacting composite fermions (as originally proposed 
by Jain g]) which can be constructed by attaching an 
even number of flux quanta to electrons under a Chern- 
Simons transformation. In this simple picture, the frac- 
tional quantum Hall effect can be mapped into the inte- 
ger quantum Hall effect for the composite fermion sys- 
tem subject to an effective mag netic field g. At filling 
factor Vj = i, although the effective magnetic field B* 
vanishes, composite fermions are subject to the random 
fluctuations of the gauge field, induced by the ordinary 
impurities Thus, it is important to study the lo- 

calization properties of non-interacting charged particles 
in a random magnetic field to understand the half-filling 
FQH system. The problem of charged particles moving 
in a random magnetic field is also relevant to the theo- 
retical studies of high T c models where the gauge filed 
fluctuations play an important role 

According to the conventional scaling theory of local- 



ization, the random flux system belongs to the unitary 
ensemble, which is described by a nonlinear sigma model 
with unitary symmetry. Since there is no net magnetic 
field, the topological term of the uniform magnetic field 
case is absent. Perturbative renormalization group calcu- 
lations show that all states are localized JhJ . However, it 
has been argued recently by Zhang and Arovas jll| that 
although the constant topological term is absent, there is 
a term describing the long ranged interaction between the 
topological densities, and they conjectured that this new 
term could lead to a phase transition from localized to 
extended states. There have been a number of conflicting 
numerical investigations on the localization properties of 
the random magnetic field system. The conclusions in 
these studies range from all states localized jl2| to ex- 
tended states around band center |@,^3|. In Ref. p4[ , 
the authors found evidence for a mobility edge, but their 
system is neither large enough to see good scaling nor 
close enough to the critical regime to obtain a conclusive 
critical exponent. 

In this paper, we systematically investigate the local- 
ization properties of a disordered two-dimensional elec- 
tron system in the presence of a random magnetic field. 
The localization length is calculated using a transfer ma- 
trix technique and finite-size scaling analysis. An impor- 
tant strength of our calculation is using system widths 
(upto 128) which are substantially (by a factor of four) 
larger than those |l2}-|l4|] existing in the literature. We 
find the following results: (1) A mobility edge E c is ob- 
served and the localization length £ diverges when the 
Fermi energy approaches the critical energy; (2) we find 
that the critical energy E c shifts with increasing disor- 
der strength; (3) E c shifts with changing randomness in 
the magnetic field; (4) the critical exponent (y w 4.5) re- 
mains unchanged while varying the disorder strength and 
the randomness of the magnetic field, indicating univer- 
sality in the metal-insulator transition; (5) the mobility 
edge survives in the presence of weak but non-zero aver- 
age random magnetic field. 

We model our two-dimensional system in a very long 
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strip geometry with a finite width (M) square lattice 
with nearest neighbor hopping. The disorder potential 
is modeled by the on-site white-noise potential Vi m (i 
denotes the column index, to denotes the chain index) 
ranging from — W/2 to W/2. Random magnetic field is 
introduced by varying the flux in each lattice plaquette 
uniformly between —<j) r /2 and <fi r /2 (in this case the av- 
erage field is zero, we also discuss the situation of weak 
but non-zero average random magnetic field in the later 
part of this paper) . The Hamiltonian of this system can 
be written as: 



A I 



(1) 



Ti = Vi 7n \im >< im\ 

i rn—1 

<im;jn> 

where < im; jn > indicates nearest neighbors on the lat- 
tice. The amplitude of the hopping term is chosen as 
the unit of the energy. A specific gauge is chosen so that 
the inter-column hopping does not carry a complex phase 
factor (i.e. ijm;i+i,m = — !)■ The only effect of random 
magnetic field shows up on the phase factor of the intra- 
column (inter-chain) hopping term. If the random flux 
in a plaquette cornered by (im), (i + 1, to), (i + 1, m + 1) 
and (i, m + 1) is 4>i m , then 
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(2) 



where <p — hc/e is the magnetic flux quantum. For a 
specific energy E, a transfer matrix Ti can be easily set 
up mapping the wavefunction amplitudes at column i — 1 
and i to those at column i + 1, i.e. 



\m(E) 
M 



= .f 



M 



(4) 



where f(x) cx 1/x in thermodynamic limit M — > oo while 
approaching a constant(~ 1) at the mobility edge where 
the thermodynamic localization length diverges. Numer- 
ically we shift the data in Fig.[l](a) onto a smooth function 
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where Hi is the Hamiltonian for the ith column, / is a 
M x M unit matrix. Using a standard iteration algo- 
rithm ||, we can calculate the Lyapunov exponents for 
the transfer matrix Ti. The localization length Xm (E) for 
energy E at finite width M is then given by the inverse of 
the smallest Lyapunov exponent. In our numerical calcu- 
lation, we choose the sample length to be over 10 4 so that 
the self-averaging effect automatically takes care of the 
ensemble statistical fluctuations. A sample of our calcu- 
lated finite width localization length for various energies 
is shown in Fig.|l|(a). 

We use the standard one-parameter finite-size scal- 
ing analysis || to obtain the thermodynamic localiza- 
tion length £. According to the one-parameter scaling 
theory, the renormalized finite-size localization length 
Am /M can be expressed in terms of a universal function 
of M/£, i.e. 



FIG. 1. (a) Renormalized finite-size localization length 
(Xm/M) for different sample width M in a random magnetic 
&e\d((p r — 1.0) without on-site disorder (W = 0); (b) Scal- 
ing function and localization length in the thermodynamic 
limit(inset) with v = 4.52 ± 0.08 and E c = -3.00. Here dif- 
ferent symbols represent different energies. 

with a least square fit. Note that we have to select data 
for a large enough sample in the scaling analysis so as 
to avoid severe finite-size effect. In our calculation, we 
choose the data for sample width greater than 16 (inclu- 
sive). The thermodynamic localization length is given 
by the amount of shifts on a log- log plot. A sample of 
the scaling function and corresponding thermodynamic 
localization length is shown in Figj^(b). Because of the 
symmetry in the problem, we only study the branch with 
negative energy E < 0. 

We first study the case with random magnetic field 
characterized by fluctuation amplitude cf> r — 1.0. We 
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find that if the Fermi energy is below E c — —3.0, the 
finite-size localization length is well converged and al- 
ways smaller than the sample width indicating that all 
states are localized below —3.0. On the contrary, for the 
electronic states with energy higher than —3.0 the in- 
verse of the Lyapunov exponent is always larger than the 
sample width which is the feature of extended states. 



FIG. 2. Scaling functions and thermodynamic localization 
length (inset) in in random magnetic fields with (b) or with- 
out (a) on-site disorder. The corresponding critical exponents 
and critical energies are: (a) v = 4.98 ± 0.10, E c = —3.13; (b) 
v = 4.86 ± 0.02, E c = -2.73. 

We also find that the thermodynamic localization length 
£ diverges while approaching E c — —3.0 indicating the 
existence of a mobility edge around —3.0. Our best fit 
analysis indeed gives a critical energy E c — —3.00 with a 
critical exponent v — 4.52 appearing in £ = \E — E c \ ~ v 
for E < E c (the critical exponent is given by the slope of 
the straight line in the inset log-log plot of Fig.[j](b)). 

We expect that the mobility edge should shift to lower 
energy if the magnetic field is less random so that ex- 
tended states are more favorable. This is exactly what 
we observe in our calculation. In Fig.|](a) we show the 
typical scaling function and localization length for a less 
random magnetic field with <p r = 0.9. In Table [| we 



present the critical exponents and critical energies for 
different randomness in the magnetic field. We find that 
the critical energy increases almost linearly with the ran- 
dom flux amplitude which is proportional to the energy 
fluctuation created by the random magnetic field. 

We now discuss the effect due to correlation between 
random magnetic field and disorder potential. We con- 
sider two types of disorder potential: (i) Independent 
model where the distribution of disorder potential is com- 
pletely independent of the random magnetic field. This 
model is relevant for the case with random distributed 
non-magnetic impurities in the sample, (ii) Correlated 
model in which the strength of disorder is associated with 
the local random magnetic field (numerically we select 
each on-site disorder so that it is proportional to the ran- 
dom flux in the neighboring plaquette) . This model is rel- 
evant for the case with random distributed magnetic im- 
purities, for example, disorder-pinned random flux lines 
in the sample. 

In the independent model, since non-magnetic ran- 
dom disorder potential tends to localize all the electronic 
states, we expect that the mobility edge in a random 
magnetic field should shift to higher energy with increas- 
ing strength of disorder potential because the localized 
states are more favorable in this situation. In Fig.^(b) we 
show the typical scaling function and localization length 
for a strongly disordered sample with W — 2.0 in a ran- 
dom magnetic field with <p r = 1.0. As presented in Table 
[jj the critical energy moves to higher energy as disor- 
der strength W increases for fixed random magnetic field 
(<j) r = 1.0) just as we expect. Presented in the last row 
of Table | is the data for a correlated disorder model, the 
behavior of the mobility edge is quite similar to that for 
the independent disorder model. 

It is interesting to notice that even though the mobility 
edge is shifting for various random magnetic field and dis- 
order strength and correlation, the critical exponent for 
the metal-insulator transition is more or less unchanged 
which indicates the universality of this critical transition. 
Our calculated exponent (y ks 4.5) is quite different 

TABLE I. The critical exponent(^) and critical energy(_E c ) 
for different representations of random magnetic field (<f> r ) and 
on-site disorder (W). The last row is for correlated disorder 
model while the other data are for the independent disorder 
model. 



W 



E c 



0.7 
0.8 
0.9 
1.0 
1.0 
1.0 
1.0 
0.8 



0.0 
0.0 
0.0 
0.0 
1.0 
2.0 
3.0 
2.0 



4.53 ±0.11 
4.60 ±0.06 
4.98 ±0.10 
4.52 ±0.08 
4.79 ±0.17 
4.86 ±0.02 
5.29 ±0.60 
4.45 ±0.10 



-3.40 
-3.28 
-3.13 
-3.00 
-2.85 
-2.73 
-2.13 
-3.35 
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from the value obtained in Ref. [Q. However, our system 
is much larger (4 times larger) than the system studied 
in Ref. [Q, our scaling is much better and our data are 
closer to the critical regime, thus our calculated exponent 
better represents the true critical exponent. 
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matrix Landauer formula type approach, finding results 
completely consistent with those presented here (i.e. an 
insulator for E < E c and a metal for E > E c ). Our ef- 
forts to construct a beta-function f3(g) have not, however, 
yielded unambiguous results because of inherent fluctua- 
tions in the calculated conductance. Systems much larger 
than those used here (M >> 128), which are not accessi- 
ble with the currently available computers, will be needed 
to obtain an unambiguous (3(g) - in that sense, we have 
a disagreement here with the conclusion of Ref. |l2| and 
|1- 
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by the U.S. National Science Foundation. 



FIG. 3. Thermodynamic localization length in the pres- 
ence of a weak but non-zero average random magnetic field. 

In Fig.^ we present the localization length for a two- 
dimensional system subject to a random magnetic field 
when the average field is weak but non-zero (i.e. < <j) >^ 
0). Clearly there exists a mobility edge at Ec ks —3.0 for 
< cj) >— 0.01. It is consistent with the argument that in 
the vicinity of < B* >= 0, the composite fermion system 
behaves as a Fermi liquid. 

As mentioned earlier there is currently considerable 
disagreement in the literature about the nature of two- 
dimensional one-electron eigenstates in a random flux en- 
vironment. An earlier finite size scaling analysis 12 con- 
cluded that all states are localized, but the localization 
lengths are exponentially large near the band center. Our 
largest system widths (=128) are four times larger than 
those (=32) used in Ref. |b| and our results are con- 
sistent with the existence of mobility edges separating 
localized and delocalized states with a localization expo- 
nent v w 4.5. We emphasize, however, that no finite- 
size scaling analysis can distinguish between extended 
states and states with extremely large (orders of magni- 
tude larger than the system length in Ref. (l2)) localiza- 
tion lengths (ljjl^l . We contend that experiments in real 
samples cannot distinguish between these two scenarios 
either i.e. the situation with extremely large localiza- 
tion lengths would behave very much like a weakly dirty 
metal. The important point to realize is that the finite- 
size scaling evidence in favor of a derealization transition 
presented in this paper is as strong as it is in the corre- 
sponding quantum Hall plateau transitions || where the 
existence of extended states is not in doubt. We have also 
calculated the conductance g as a function of the system 
length for various (Fermi) energies using a direct transfer 
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